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Abstract. Similar directions are obtained for the local interstellar magnetic field (ISMF) by com- 
paring diverse data and models that sample five orders of magnetic in spatial scales. These data 
include the ribbon of energetic neutral atoms discovered by the Interstellar Boundary Explorer, he- 
liosphere models, the linear polarization of light from nearby stars, the Loop I ISMF, and pulsars 
that are within 100-300 pc. Together these data suggest that the local ISMF direction is correlated 
over scales of ~ 100 pc, such as would be expected for the interarm region of the galaxy. 
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1. INTRODUCTION 

The properties of magnetic turbulence in interarm versus arm regions are distinctly 
different, as shown by Faraday rotation measure (RM) structure functions determined 
separately for the two regimes. The interstellar magnetic field (ISMF) is coherent over 
100 pc scales in interarm regions, and random over spatial scales of less than 10 pc in 
arm regions [1]. Which of these scenarios dominates for the ISMF in the solar vicinity? 
The heliosphere, and its ribbon of energetic neutral atoms (ENAs) discovered by the 
Interstellar Boundary Explorer (IBEX) spacecraft [2], probe the direction and strength 
of the very local ISMF in a single spatial location. Over intermediate scales, the position 
angles of starlight polarized by magnetically aligned interstellar dust grains trace the 
ISMF direction. Pulsar RM data give the field direction for spatial scales of 100-300 
AU. As discussed in Frisch et al. [3], together these data and cosmic ray anisotropics 
suggest the ISMF near the Sun is similar to expectations for interarm regions. 

2. IBEX RIBBON, HELIOSPHERE, AND LOCAL ISMF 



IBEX measures ENAs formed by charge-exchange (CEX) between heliospheric ions 
and neutral gas from the surrounding ISM, which flows through and around the helio- 
sphere at a velocity of ~ 26 km s _1 from the direction £,b~ 4°, 15°. Among the sources 
of the ions creating the 0.01-6 keV ENAs detected by IBEX are the post-shock solar 
wind, pickup ions, and non-thermal energetic ions in the inner heliosheath, secondary 
ions created by CEX between interstellar protons and ENAs escaping beyond the he- 
liopause, and possibly pickup ions inside of the solar wind termination shock. The first 
ENA full skymap, collected early 2009, revealed an unexpected "ribbon" of ENAs that 
traces heliosphere asymmetries created by the angle between the ISMF and velocity of 




FIGURE 1. (A) Contours of constant angular offsets from the pole of the IBEX ribbon at A=221°, 
/3=39°, shown plotted against the IBEX 1 keV ENA fluxes. The black line shows the locus of points 
perpendicular to the ISMF direction, e.g. BR=0, based on the ISMF displayed in (C) from the heliosphere 
model of [1 1]. B is the magnetic field vector and R is the radial sightline. Plot (B) shows the correlation 
between global 1 keV ENA fluxes and B R. Plot (A) is in ecliptic coordinates, and centered on the 
heliosphere nose. The locations of the two Voyager spacecraft are identified. The figure is from [5]. 



inflowing neutral interstellar gas [Fig. I, 2, 4, 5]. The ribbon is visible in the energy 
range ~ 0.2 — 6 keV and forms a nearly complete circle in the sky, that is offset by 
~ 75° from the ribbon arc center [6]. Comparisons between heliosphere models (see 
below) and the ENA ribbon revealed that the ribbon is seen towards sightlines that are 
perpendicular to the ISMF draping over the heliosphere [5]. Seven possible formation 
mechanisms for the ribbon have been discussed [7, 6, 2, 8]. In order to link the ribbon di- 
rectly to the ISMF draping over the heliosphere, the mechanism based on CEX between 
interstellar H° and energetic ions beyond the heliopause is assumed to be correct [9, 10]. 
The ribbon location is stable over the first several skymaps [7]. The ISMF direction at 
the heliosphere is then given by the center of the ribbon arc, located at £=33°, b=55° (see 
Table 1 for ecliptic coordinates). 

The best match between the predicted and observed ribbon locations occurs for 
ENAs produced roughly 100 AU beyond the heliopause [12]. The surrounding cloud 
is ~ 22% ionized [n(p) ~ 0.07 cm~ 3 , 13]. The ribbon formation mechanism requires 
that ENAs escaping the heliosphere CEX with interstellar protons piled up against 
the heliopause, creating a secondary pickup ion population in the outer heliosheath. 
The secondary pickup ion must CEX quickly with an interstellar H° atom (before 
magnetic turbulence disrupts the ring-beam distribution), to create an ENA detected 
preferentially in directions perpendicular to the ISMF The relative roles of large- and 
small-scale magnetic turbulence in the outer heliosheath are not presently understood, 
so this mechanism is still under study [14, 15]. A recent model [12] reproduces the 
ribbon for an ISMF strength of < 3 /lG, with the pole directed from £=33°, b=54° 
(or A,/3 = 222°, 39°). A local ISMF strength of ~ 3 /lG is found from photoionizaton 
models of the local interstellar cloud (LIC) around the heliosphere, and the assumption 



of pressure equilibrium between the ISMF and thermal gas [13]. 

MHD heliosphere models are constructed to reproduce the heliospere asymmetries 
induced by the ~ 45° angle between the ISMF direction and gas flow direction [11, 16, 
17, 18]. The ISMF direction in these models is generally directed downwards through 
the ecliptic plane in order to reproduce asymmetries such as the ~ 5° offset between 
the central directions of the flows of interstellar H and He into the heliosphere, and the 
difference between the solar wind termination shock distances encountered by Voyager 
1 (94 AU) versus Voyager 2 (84 AU). 

If the IBEX ribbon is formed beyond the heliopause through CEX with secondary 
pickup ions, then the ribbon is an extraordinary diagnostic of the local ISMF and 
neutral densities. Models show that variations in interstellar properties of ~ 15% lead to 
pronounced differences in the location and width of the ribbon because of the interplay 
between magnetic, plasma and neutral pressures [19, 12]. 

Two aspects of the ribbon data are notable: (1) The imprint of the solar wind on ENA 
fluxes and the ribbon is apparent in flux decreases of ~ 15% between the first and second 
skymaps, spaced by six months [7] . The travel time for a 1 keV solar wind proton to pass 
the Earth, the solar wind termination shock, and then return to IBEX as an ENA, is over 
two years. This time lag suggests that the decrease in 0.7-1.1 keV fluxes in the second 
map is from the unusually low solar wind pressure of the recent solar cycle minimum. (2) 
The energy distribution of global ENAs shows an overall power law spectra, E~ K , with 
K ranging between 1.28 and 2.04 for different regions of the sky and spectral intervals 
[4]. The power-law distribution of ribbon ENAs becomes softer above a "knee" in the 
spectrum between 1-4 keV, which occurs at higher energies for higher latitudes [6]. 

3. LOCAL ISMF DIRECTION FROM POLARIZED STARLIGHT 

The galactic ISMF was first mapped in the 1940's using polarized starlight, which 
showed patterns of linear polarization related to dust extinction. In the diffuse ISM, dust 
grain alignment occurs when magnetic torques align asymmetric rotating charged dust 
grains, producing a birefringent interstellar medium with lowest opacities parallel to the 
ISMF direction. This gives linear polarization that is parallel to, and traces, the ISMF 
direction [20]. For sightlines perpendicular to the ISMF direction in a uniform medium, 
polarization will increase with distance. Polarizations do not provide the strength or 
polarity of the ISMF, and polarization is reduced by foreshortening towards the poles of 
the field. 

During the early 1970's, Tinbergen [22] conducted a sensitive (for that time) study 
of the interstellar polarizations of stars within ~ 40 pc of the Sun, and found a patch of 
dust in the south and close to the Sun. Based on interstellar absorption lines towards a 
Cen at 1.3 pc, and polarization towards 36 Oph, the magnetically aligned dust must 
be in the "G-cloud" within 1.3 pc [23]. Recent polarization data confirm that local 
interstellar polarizations are significantly weaker in the northern hemisphere than the 
southern hemisphere. 

We have analyzed available interstellar polarization data, including six recent sets of 
sensitive data, to determine the best fitting ISMF direction for a 180° diameter region 
that is centered on the heliosphere nose and coinciding roughly with the G-cloud ([21, 3]. 



mean(sin(PA j )) i and B grld = B i=besl =263°, 37° 




FIGURE 2. Best best-fitting ISMF direction from polarization data is determined by minimizing the 
function < sin (PA;) >. PA; is the position angle of linear polarizations in the i th possible ISMF direction. 
The mean is evaluated for data towards stars within 40 pc and 90° of the heliosphere nose (where the G- 
cloud is seen, see [21]). Gray dots show the best fitting ISMF, which is directed towards A=263°, /3=37°, 
with uncertainties of ~ ±35° (corresponding to £,b=38°, 23°). A new study using an extended polarization 
data set returns similar results for the local ISMF direction, but the uncertainties are larger [3] . Polarization 
vectors are plotted with either black or white bars. The figure is from [21]. 



Although we utilize polarization towards stars out to 40 pc, most of the polarization 
should be formed in the region with the highest concentrations of gas within ~ 15 pc 
[23, 3]. These data indicate a local ISMF direction with a pole towards £,b~ 38°, 23°, 
with uncertainties of ~ ±35° (Fig. 2, Table 1, [21, 3]). This direction is centered 
approximately 35° from the center of the IBEX ribbon arc. The small difference between 
the two directions may indicate turbulence in the local ISMF, for example between 
the ISMF directions of the G-cloud (within 1.3 pc) and the LIC in which the Sun is 
embedded. In several regions, the flow of interstellar gas past the Sun shows merging or 
colliding clouds [24], such as the two clouds in front of Sirius [25], which could drive 
local magnetic turbulence. 



TABLE 1. Interstellar Magnetic Field Direction, 150 AU - 150 pc 



Scales 


ISMF Source 


Angle 


A,j3 




- 150 AU 


IBEX ribbon center [4] 


0° 


221°, 39° 


33°, 55° 


- 150 AU 


Heliosphere model [12]* 


3° 


224°, 41° 


36°, 53° 


< 40 pc 


Optical polarization [21, 3] 


32° 


263°, 37° 


38°, 23° (±35°) 


80 pc 


SI Magnetic Bubble [26] 


46^1 


291°, 64° 


71°, 18°(±48°) 


150 pc 


Pulsar RM data [27] f 


23° 


232°, 18° 


5°, 42° 



* The ISMF polarity is directed downwards through the ecliptic plane 
' The ISMF polarity is directed upwards through the ecliptic plane 



4. THE LOOP I SUPERBUBBLE AND LOCAL ISMF 



Loop I is an evolved but reheated superbubble that dominates the northern hemisphere 
because of it's large angular radius. It appears as a distinct imprint on the ISMF within 
100 pc that is detected in polarized starlight, polarized radio continuum, Faraday rota- 
tion, and Zeeman splitting in the HI filaments defining the shell circumference [28]. All 
existing spherically symmetric models of the radio continuum Loop I place the Sun in 
the superbubble rim. A recent study of the intensity of 1.4 GHz emission in this region 
derived two subshells of Loop I containing enhanced synchrotron emissivity from swept 
up magnetic fields, and placing the Sun in the rim of the SI subshell [26, 29]. The direc- 
tion of the bulk flow of the "cluster of local interstellar clouds" (CLIC) past the Sun has 
an upwind direction of £,b=335° -7° ([25], after correcting for the solar velocity through 
the local standard of rest, LSR), near the center of the S 1 shell, 

If the CLIC is part of the nearside of an evolved superbubble shell associated with 
Loop I [28, 23], then both the polarity and direction of the local ISMF could be ordered 
over the ~ 78 pc radius of the S 1 subshell. The angle between the local ISMF direction 
and the LSR CLIC velocity is 68°-78°, depending on whether the comparison is with 
the polarization axis or ribbon arc center. This large angle suggests that the very local 
ISMF is roughly parallel to the local surface of the expanding Loop I shell. 

5. ISMF FROM PULSARS IN THE LOCAL BUBBLE 

The ratio of pulsar Faraday rotation measures (RM) and dispersion measures provides an 
electron-density weighted measurement of the ISMF component parallel to the sightline, 
and the field polarity. Salvati (2010) performed a best- fit to the rotation and dispersion 
measures of four pulsars within 160-300 pc, and in the very low density third galactic 
quadrant corresponding to the interior of the Local Bubble. He found a best-fit direction 
of £~ 5°, b~ 42°, and a field strength of |fi|=3.3 \lQ (Table 1). This direction is within 
23° of the center of the IBEX ribbon arc, and 33° of the best-fitting ISMF direction. The 
polarity of this ISMF is directed towards the pole of A,/3=232°,18°, indicating that it is 
directed upwards through the plane of the ecliptic. 1 Wolleben [26] finds a similar RM 
polarity in this direction, although variations over small angular scales are seen. 

6. GALACTIC COSMIC RAYS AND THE LOCAL ISMF 

An unusual indicator of the ISMF affecting the heliosphere is provided by GeV-TeV 
galactic cosmic rays (GCR). Cosmic rays at sub-TeV energies have gyroradii less than 
several hundred AU in a 1 /iG magnetic field, which is comparable to heliosphere 
dimensions. GCRs with energies < 10 3 GeV show pronounced spatial asymmetries 
[30, 31, 32], including a sidereal distribution with a southern GCR deficit, and a broad 
excess emission region attributed to the heliotail. The broad tail-in excess at 500-1000 



Note that the rotation measure is positive when the ISMF is directed towards the observer, by definition. 



GeV is centered at ecliptic coordinates of A,j8 of 71°,-42° to 90°,-53°, and is skewed 
with respect to ecliptic latitude and covers the direction opposite to the ISMF defined by 
the ribbon arc (A,/3=41°,-39°, Table 1). This tail-in excess also coincides with the axis 
of the ISMF determined from the polarization data, A,/3=83°,-37°. The minimum in tail 
ENA fluxes is 44° west of the downwind gas-flow direction [6], and it overlaps the excess 
of sub-TeV particles. Evidently the GCR excess could form either in deep tail regions or 
indicate GCRs streaming along the local ISMF or SI shell. Lazarian and Desiati (2010) 
have modeled the tail-in sub-TeV excess as due to the stochastic acceleration of particles 
in the time-variable magnetic field of the heliotail. 

7. SUMMARY 

The summary of data in the previous sections shows that the direction of the ISMF ap- 
pears to be relatively constant over scales of five orders of magnitude, with warping 
(large-scale turbulence) or variations in the field direction typically less than 30°-40° 
(Table 1). The agreement between ISMF direction from pulsars in the Local Bubble, 
starlight polarizations, the IBEX ribbon, and heliosphere models, is remarkable consid- 
ering that the spatial scales of these estimates differ by five orders of magnitude. The 
puzzle is that opposite polarities are inferred for the ISMF shaping the heliosphere and 
the pulsar RM data, downwards through the ecliptic plane in the first case and upwards 
through the ecliptic plane in the second case. It could be pure coincidence that the Local 
Bubble ISMF direction from the pulsar data and from heliosphere measurements and 
models are similar. However the galactic magnetic field is known to be ordered over 
kiloparsec spatial scales in low density interarm regions of the Galaxy such as around 
the Sun [33]. Perhaps it is not so surprising to see uniformity of the magnetic field di- 
rection over local scales of 100 pc. If the ISMF traced by the IBEX ribbon, polarization 
data, Loop I shell, and pulsar data are related, it would mean that Loop I is a highly 
asymmetric evolved superbubble, with the Sun in the segment most distant from the 
superbubble source and most closely linked to the interarm field direction. 
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